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FOREWORD 

The  work  described  in  this  report  was  carried  out  under 
Contract  NAS  1-3778  with  the  National  Aeronautics  and  Space 
Administration. 

This  is  a  revision  of  the  original  report,  which  was  issued  in 
August  1964. 
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SUMMARY 


(Recent  investigations  by  Stein  and  by  Fischer  on  the  influence  of  edge  conditions  on 
the  critical^load'of  cylindrical  shells  are  here  extended  to  cover  six  additional  com¬ 
binations  of  boundary  conditions.  The  results  show  that  drastic  reductions  of  the 
critical  load  for  cylinders  with  lateral  support  of  the  edges  are  obtained  only  if  the 
edges  are  free  in  the  tangential  direction.  For  other  boundary  conditions,  this  reduc¬ 
tion  is  never  more  than  about  20  percent.  Consequently,  the  results  of  this  investiga¬ 
tion  alone  cannot  explain  the  well-known  discrepancy  between  theory  and  test  data. 
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NOTATION 


integration  constants;  see  Eqs.  (19)  and  (20) 

Et3/[  12(1  -  v2)] 

Young's  modulus 

stress  function;  see  Eq.  (22) 

value  of  F  at  x  =  xm 

m 

shell  length 

total  number  of  discrete  points 
compressive  axial  load  per  unit  width 
N/(2y) 

axial,  circumferential,  and  shear  forces  per  unit  width 
lateral  displacement;  see  Eq.  (22) 


stress  function 


parts  of  stress  function,  prebuckling  and  incremental 
distance  between  neighboring  discrete  points 


number  of  waves  in  circumferential  direction 


internal  pressure 


p  (pr)/(Et) 


r 

t 


U,  V,  w 


w 


O’ 


w 


1 


X,  y 


X 


in 


m 


a 


a 


ex’  V  yxy 


shell  radius 
shell  thickness 

axial,  circumferential,  and  lateral  displacements 

parts  of  lateral  displacement,  prebuckling  and  incremental 

axial  and  circumferential  coordinates 

value  of  axial  coordinate  at  discrete  point 

see  Eq.  (28) 

(a1L)/(2r) 

(a2L)/(2r) 

axial,  circumferential,  and  shear  strain 
(t/ r)/[  12(1  -  ^2)]1/2 


ctCR 

ctCL 

4 


Poisson's  ratio 
critical  axial  stress 

classical  value  of 

4,4  4,22  4,4 

a  /ax  +  2a  /ax  ay  +  a  /ay 


When  subscripts  x  and  y  follow  a  comma,  they  indicate  partial  differentiation  of  the 
principal  variable  with  respect  to  x  or  y  .  Primes  indicate  total  derivatives  with 
respect  to  x  . 


x 


Section  1 
INTRODUCTION 


Serious  disagreement  between  the  results  of  classical  theory  and  experiments  on  the 
buckling  of  axially  compressed  cylindrical  shells  has  long  been  known  to  exist.  Despite 
considerable  effort,  this  problem  is  still  not  completely  understood.  New  light  was 
thrown  on  the  matter  in  the  recent  investigation  by  Stein  reported  in  Refs.  1  and  2. 

In  that  analysis,  the  deformations  and  stresses  induced  by  edge  support  were  for  the 
first  time  rigorously  taken  into  account.  It  was  shown  that  previous  inconsistent 
assumptions  in  the  classical  theory  with  regard  to  boundary  conditions  in  the  prebuckling 
state  can  significantly  affect  the  theoretical  buckling  load. 

The  analysis  of  Ref.  1  is  limited  to  one  specific  set  of  simple-support  boundary  con¬ 
ditions.  Fischer,  in  Ref.  3,  later  considered  the  same  problem  but  used  different 
in-plane  boundary  conditions.  The  results  of  the  two  analyses  differ  significantly,  and 
it  appeared  desirable  that  a  more  complete  investigation  be  undertaken  to  establish  the 
effect  of  edge  conditions  on  the  classical  buckling  load  for  axially  compressed  cylindrical 
shells.  Such  an  analysis  is  presented  here  and  includes  cylinders  with  clamped  as 
well  as  simply  supported  edges.  Eight  different  sets  of  boundary  conditions  are  con¬ 
sidered,  and  the  approach  to  the  problem  is  largely  the  same  as  in  Ref.  1.  The  pre¬ 
buckling  displacements  are  obtained  through  rigorous  solution  of  the  applicable  dif¬ 
ferential  equations.  As  incremental  displacements  are  assumed  to  be  infinitesimal, 
the  buckling  equations  are  linear  but  have  variable  coefficients.  Solutions  to  these  are 
obtained  by  use  of  a  finite  difference  technique. 

It  was  pointed  out  by  Koiter  (Ref.  4)  that  the  classical  solution  of  the  buckling  problem 
for  axially  compressed  cylinders  is  based  on  the  boundary  conditions  used  by  Fischer, 
and  that  Stein's  condition  of  zero  tangential  edge  force  may  lead  to  a  reduced  critical 
load  even  if  a  membrane  solution  is  used  for  prebuckling  displacements.  This  was 


shown  to  be  the  case  by  Ohira  (Ref.  5)  and  by  Hoff  and  Rehfield  (Ref.  6).  Therefore, 
the  present  analysis  includes,  for  comparison,  solutions  for  the  case  in  which  effects 
of  edge  restrain  on  prebuckling  configurations  have  been  omitted.  This  case  will  be 
referred  to  as  the  membrane  prebuckling  solution. 
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Section  2 

BASIC  EQUATIONS 


As  in  Ref.  1,  Donnell -type  equations  are  used  in  the  analysis.  In  their  nonlinear 
form  the  three  equilibrium  equations  are 


N  +  N  =  0 
x,  x  xy,y 


N  +  N  =  0 
y.y  xy,x 


(1) 


DV  w  +  N/r-(Nw  +Nw  +2N  w  )  -  p  =  0 
y  '  x  ,xx  y  ,yy  xy  ,xy' 

The  corresponding  relations  between  stresses,  strains,  and  displacements  are 


Nx  =  [Et/(1  -  v  )]  (ex  +  vey) 


Ny  =  [Et/(1  -  v)]  ( ey  +  vej 


(2) 


N  =  [Et/2( 1  +  i>)]  y 
xy  v  xy 


e  =  u  +  1/2  w 
X  ,x  ,x 


=  v  „  +  w/r  +  1/2  w 

y  »y  *y 


(3) 


y  =  u  +  v  +ww 
xy  ,y  ,  x  ,x  ,y 
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The  first  two  of  the  equilibrium  equations  are  identically  satisfied  if  a  stress  function, 
f  ,  is  introduced  such  that 


N  =  f  ;  N  =  f  ;  N  =  -  f 
x  ,yy  y  ,xx  xy  ,xy 


By  use  of  Eqs.  (4),  the  third  equilibrium  equation  can  now  be  written 


(4) 


DV4w 


w 

,yy  »xx 


f  w 

,  XX  , 


Through  elimination  of  u  and  v  in  Eqs.  (2)  and  (3),  a  compatibility  equation  can  be 
derived: 


( l/Et )  V 


+  w  w 
,  xx 


Equations  (5)  and  (6)  are  two  simultaneous  differential  equations  in  the  two  unknown 
variables  f  and  w  .  Together  with  an  adequate  set  of  boundary  conditions,  they 
govern  the  behavior  of  the  shell.  A  total  of  eight  boundary  conditions  -  four  at  each 
edge  of  the  cylinder  —  is  required.  As  the  displacements  are  expected  to  be  sym¬ 
metrical  about  the  midpoint  of  the  shell,  it  is  also  possible  to  use  a  set  of  four  sym¬ 
metry  conditions  together  with  four  conditions  at  the  edge  of  the  shell. 


The  axial  coordinate  is  assumed  to  be  zero  at  the  midpoint  and  consequently 


at  x  =  0  :  w  =  0 

,  x 


w 

,  XXX 


0 


f  =  0 

,  XXX 


(7) 
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The  shell  is  assumed  to  be  supported  in  the  radial  direction  at  the  edge.  Therefore, 

at  x  =  L/2:  w  =  0  (8) 

Three  additional  conditions  are  needed  and  will  be  obtained  by  selection  of  one  of  each 
of  the  following  three  pairs: 

at  x  =  L/2  :  w  v  =  0  or  w  =  0 

u  —  0  or  N  =  0  (9) 

x 

v  =  0  or  N  =0 

xy 

Evidently,  eight  different  combinations  are  possible  and  these,  in  combination  with 
Eqs.  (7)  and  (8),  form  eight  complete  sets  of  boundary  conditions.  All  of  these  sets 
will  be  considered  in  the  analysis. 

For  shorter  writing,  the  different  combinations  will  be  referred  to  as  follows: 


Case  SI 

when 

w 

,  XX 

=  0  , 

U  =  0  , 

V 

=  0 

S2 

w 

,  XX 

=  0  , 

o 

II 

V 

=  0 

S3 

w 

,XX 

=  0  , 

u  =  0  , 

N 

xy 

=  0 

S4 

w 

,  XX 

=  0  , 

N  =  0  , 

X 

N 

xy 

=  0 

Cl 

w 

,x 

=  0  , 

u  =  0  , 

V 

=  0 

C2 

w 

,X 

=  o  , 

O 

II 

V 

=  0 

C3 

w 

,X 

=  0  , 

u  =  0  , 

N 

xy 

=  0 

C4 

w 

,X 

=  0  , 

N  =  0  , 

X 

N 

xy 

=  0 
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Section  3 

PRE BUCKLING  DISPLACEMENTS 


Displacements  and  stresses  will  be  composed  such  that 


w  =  Wq  + 


f  ■  *0  *  fl 


where  subscript  0  identifies  the  conditions  at  impending  buckling  and  hence  subscript 
1  identifies  incremental  quantities. 

At  impending  buckling  the  deformations  are  axially  symmetrical  and  therefore  the 
equation  of  equilibrium  can  be  written 


Dwo,xxxx  +  Wr  '  f0,yyw0,xx  -  0  ’  0 


If  the  applied  compressive  axial  load  is  denoted  by  N 


f„  =  -  N 

o.yy 


f0,xx  =  Etw0/r  -  ,'N 


Substitution  of  Eqs.  (12)  into  Eqs.  (11)  yields 


Dw0,xxxx  +  Nw0,xx  +  Etw0/r  -  <1/N/r  +  P)  -  0 
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The  solution  to  this  equation  can  be  written  in  the  form 


wQ  =  r(2j/yN  +  p)(l  +  A1  sin  a1x  sinh  a2x  +  cos  axx  cosh  a2x)  (14) 


where 

al  = 

k(  1  +  N)1//2 

a2 

k(l  -  N)l/2 

k  = 

(2y  )_1/2 

r  = 

(t/r )/[l2( 1  -  v2)\  f 

N  - 

N/(2y) 

p  = 

(pr)/(Et) 

(15) 

The  integration  constants  are  obtained  through  substitution  of  Eq. 

(4)  into  the  boundary 

conditions. 

For  simply  supported  shells,  these  are 

at  x  =  L/2:  wQ  =  0  ,  wQ>xx  =  0  (16) 


and  for  clamped  shells 

at  x  =  L/2:  wQ  =  0  ,  wQ  x  =  0  (17) 

The  following  notations  are  introduced. 

=  (a1L)/(2r) 
a2  =  (a2  L  )/(2r ) 


(18) 


For  simply  supported  shells  the  solution  will  be  given  by 


r  — 2  l/2  _  I 

Ai  =  -  (1  -  N  )  sin  a^  sinh  a2  +  N  cos  cosh  of  2  I  / 

(1  -  N2)1//2  (cosh2  a  £  -  sin2  Ck!  ^  )  j 
A2  =  ~  [(-1  "  N2)172  cos  a1  cosh  a2  -  N  sin  a  1  sinh  a 2  ]/ 


(1  -  N2)1//2  (cosh2  a  2  -  sin2  a.jjJ 


(19) 


For  clamped  shells 


A.^  =  -  2  |(1  -  N)1//2  cos  a i  sinh  a2  -  (1  +  N)^2  sin  a^  cosh  a2j  / 

|(1  +  N)1//2  sinh  2a 2  +  (1  -  N)1//2  sin  2a jJ 
A2  =  -  2  |(  1  -  N)*^2  sin  a L  cosh  a2  +  (1  +  N)  ^  cos  a^  sinh  a2J  / 

|(1  +  N)1^2  sinh  2a2  +  (1  -  N)1//2  sin  2a.J  (20) 
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Section  4 

BUCKLING  EQUATIONS 


The  equations  needed  for  the  solution  of  the  buckling  problem  are  obtained  through 
substitution  of  Eqs.  (10)  into  Eqs.  (5)  and  (6).  From  the  equations  so  derived,  the 
prebuckling  equations  can  be  subtracted  and  higher  order  terms  in  the  infinitesimal 
incremental  displacements  omitted.  The  equilibrium  and  compatibility  equations  may 
then  be  written  in  the  form 


DVS  +  Wr  +  Nwl.xx  -  "Wl.yy  +  <rN  "  Etw0/r)wl, 


yy 


( l/Et  )V  fx-  w1)XX/r  +  w0)XXw1( 


yy 


o 

(21) 

0 


The  equations  may  be  separated  with  respect  to  the  space  variables  by  use  of  the 
substitution 


=  F(x)  sin  (ny/ r) 
wj  =  W(x)  sin  (ny/r) 


(22) 


The  conditions  of  continuity  in  the  circumferential  direction  will  be  satisfied  if  n  , 
the  number  of  circumferential  waves,  is  an  integer.  In  substitution  of  Eqs.  (22)  into 
Eqs.  (21),  primes  are  used  to  denote  differentiation  with  respect  to  x  . 

D[W""  -  2(  n/r  )2 W"  +  (n/r)4W]  +  F"/r  +  NW"  +  (n/r)2w'^F 

-  (n/r)2  (vN  -  EtwQ/r)W  =  0  (23) 

(l/Et)  [F""  -  2( n/r )2 F"  +  (n/r)4F]  -  W"/r  -  (n/r)2w|JW  =  0 


11 


These  are  linear  differential  equations  with  variable  coefficients.  Besides  the  trivial 
solution  W  =  F  =  0  ,  solutions  which  satisfy  the  boundary  conditions  exist  for  partic¬ 
ular  values  of  N  .  The  lowest  of  these  values  represents  the  critical  load  of  the  . 
cylinder. 

The  boundary  conditions  are  given  in  Eqs.  (7),  (8),  and  (9).  The  conditions  pertaining 

to  w  f ,  N  ,  or  N  are  converted  in  terms  of  W  and  F  in  an  obvious  way.  The 
’  x  xy 

conditions  for  u  and  v  may  be  expressed  in  terms  of  W  and  F  by  use  of  Eqs.  (2), 
(3),  (4),  and  (22). 

It  is  found  that  the  condition  v  =  0  can  be  replaced  by 

F"  +  r(  n/r  )2F  =  0  (24) 


and  the  condition  u  =  0  by 


F»"  -  (2  +  v)  (n/r)2F'  -  EtW'  =  0  (25) 

The  classical  formulation  of  the  problem  is  obtained  if  the  membrane  solution 

WQ  =  ( i>Nr )/(  Et )  +  (pr2)/(Et)  (26) 


is  substituted  in  Eqs.  (23). 

The  simplified  equations  so  obtained  for  the  case  of  hydrostatic  pressure  were  solved 
by  Sobel  in  Ref.  7  for  the  eight  sets  of  boundary  conditions  considered  here.  Cor¬ 
responding  solutions  presented  here  for  the  case  of  axial  compression  were  obtained 
by  use  of  Sobel's  computer  program.  These  solutions  are  compared  with  solutions 
of  Eqs.  (23). 
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Section  5 

METHOD  OF  SOLUTION 


Solutions  to  the  equations  governing  the  buckling  problem  are  obtained  through  use  of 
the  finite  difference  method.  The  equations  are  approximately  satisfied  at  a  number  of 
discrete  points  along  the  shell  length  in  the  interval  0  <  x  ^  L/2  .  In  order  that  the 
eight  boundary  conditions  may  be  defined,  four  additional  points,  two  at  each  end,  are 
located  outside  the  interval.  The  points  are  chosen  to  be  equally  spaced,  and  the  dis¬ 
tance  between  adjacent  points  is  denoted  by  h  .  The  total  number  of  points,  hence,  is 

M  =  L/(2h)  +  5 

The  axial  coordinates  of  the  points  are  given  by 

xm  =  h(m  -  3)  ;  m  =  l,2,3,...,M 

The  values  of  the  stress  function  and  the  lateral  displacement  at  the  point  with  the 

axial  coordinate  x  are  referred  to  as  Fm  and  .  At  each  of  the  points, 
m  mm 

F  and  W  are  unknown  and  hence  the  total  number  of  unknowns  is 
m  m 

2M  =  ( L/h  +  10 ) 

At  each  of  the  interior  points  (m  =  3 , 4  . . . , M  -  2)  the  compatibility  and  equilibrium 
equations  are  formulated  in  terms  of  finite  differences,  yielding  2[L/(2h)  +  1]  = 

L/h  +  2  equations.  The  boundary  conditions,  written  in  terms  of  finite  differences, 
provide  eight  equations  for  a  total  of  (  L/h  +  10 )  . 


13 


The  following  difference  approximations  are  used: 


y|  =  iyi+i  -  yi-!1'''2") 

y'i  ■  iyi+i  -  2yt  +  1/1,2 

yj"  -  u/2yM  -  y1+i +  n-i  - 1/2  ^i-z ) 7113 
yj"’  =  iy1+2  -  4yI+i + ^  -  4yi-i +  yi-2 1/1,4  (27 

In  order  that  nonzero  elements  in  the  coefficient  matrix  will  be  concentrated  around 
the  main  diagonal,  the  first  four  equations  are  the  boundary  conditions  at  midlength, 
and  the  last  four  the  boundary  conditions  at  the  edge.  Furthermore,  of  the  remaining 
equations,  the  odd  numbered  are  compatibility  equations  and  the  even  numbered  are 
equilibrium  equations.  The  unknowns  are  denoted  y_.  ,  such  that 

y2m  ^m 
y2m-l  Fm 


where 

1  <  m  <  L/(2h)  +  5  (28) 

It  follows  from  Eqs.  (27)  and  (28)  that 

Fm  =  [y2m+l  "  y2m-3]/(2h) 

Fm  -  ^y2m+l  -  2y2m-l  +  y2m-3^//^1 

Fm  =  ^ 1//2  y2m+3  ~  y2m+l  +  y2m-3  "  y2m-5  ^  ^ 

Fm  _  ^y2m+3  _  4y2m+l  +  ^y2m-l  "  4y2m-3  +  y2m-5 
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Wm  '  ly2m+2  -  y2m-2l/(2h> 

Wm  =  ly2m+2  -  2y2m  +  y2m-2  1/1,2 
Wm  =  ll/2y2m+4  ‘  y2m+2  +  y2m-2  ’  1/2  y2m-4  1/1,3 
Wm"  -  ly2m+4  ‘  4y2m+2  +  6y2m  “  4y2m-2  *  y2m-4 1/1,4 
By  use  of  Eqs.  (29),  the  boundary  conditions  at  midlength  may  be  written 


Equation  1  ; 

F3 

=  (y7  -  y3 )/2h  =  o 

2  ; 

pm 

3 

=  ( 1/2  y9  -  y7  +  y3  -  */2  y^/1*3 

3  ; 

W3 

=  (y8  -  y4)/2h  =  o 

4  ; 

W3U 

=  ( 1/2  y  10  -  y8  +  y4  -  i/2y2)/h' 

The  compatibility  conditions  provide  the  equations  numbered  5,7,9,...,  2M  5  . 

For  m  =  3,4,5,...,  M-2  , 

(l/Et)  [(y2m+3  "  4y2m+l  +  6y2m-l  ~  4y2m-3  +  y2m-5 

-  2<n/r>2  ly2m+l  -  2y2m-l  +  y2m-3)/h2  +  (n/r ,4y2m-ll 

-  <1/rl  <y2m+2  -  2y2m  *  y2m-2>/h2  ‘  (n/r)2  (W0>my2m  =  ° 

(31) 


Similarly  the  equations  numbered  4,6,8,...,  2M-4  are  obtained  from  equilibrium 
conditions. 


15 


For  m  =  3,4,5,...  ,M-2  , 


D  I  ^  y2m+4  -  4y2m+2  +  6y2m  ‘  4y2m-2  +  y2m-4)/h4  ‘  2(n/r)Z  <y2m+2  ‘  2y2m 

+  y2m+2)/h2  4  (n/r)4y2ml  +  <l/r)(y2m+l  '  2y2m-l  +  y2m-3  ^  *  N(  y2m+2 

-  2y2m  +  y2m-2)/h2  +  <n/r>2<w5>  y2m-l  -(n/r,Z  I"14  ‘  Et(w0)  1  y2m  *  0 

m  mJ 

(32) 

The  last  four  boundary  conditions  in  terms  of  finite  differences  are 


Equation  2M-3: 

3 

ii 

o 

’  y2M-4 

0 

Equation  2M-2: 

o 

II 

X 

’  *y2M-2 

y2M-6^2h  =  ° 

or 

W  =  0 

,xx 

’  (y2M-2 

2y2M-4  +  y2M-6^h  =  ° 

Equation  2M-1: 

u  =  0 

’  ^1//2y2M- 

1  "  y2M-3  +  y2M-7  "  1/<2y2M-9^h3 

(2  +  v )  ( n/r  )  (y2M_3  ~  y2M-7^2h^ 

Et(  y2M 

-2  "  y2M-6^2h^  0 

or 

N  =  0  : 

X 

;  y2M-5 

0 

Equation  2M: 

v  =  0  ; 

;  '  y2M-3 

2y2M-5  +  y2M-7)/<h2)+  1'(”/r)2y2M- 

or 

o 

n 

& 

z 

1  ^  y2M-3 

y2M-7>/2h=° 

Once  the  boundary  conditions  have  been  chosen,  Eqs.  (30)  through  (33)  define  a  linear 
homogeneous  equation  system  of  the  order  2M.  This  equation  system  can  have  nontrivial 
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solutions  only  if  the  determinant  of  the  coefficient  matrix  equals  zero.  The  lowest 
value  of  N  for  which  the  determinant  is  zero  represents  the  critical  load  of  the  cylinder. 
This  value  is  found  through  computation  of  the  determinant  for  a  series  of  values  of  the 
load  such  that  the  solution  can  be  found  graphically. 

For  the  computations,  an  IBM-7094  digital  computer  was  used.  In  programming, 
advantage  was  taken  of  the  fact  that  elements  sufficiently  far  off  the  diagonal  were 
zero.  Furthermore,  the  computer  program  provided  a  possibility  to  alternately 
reduce  the  matrix  and  compute  new  coefficients.  In  this  way  the  available  storage 
space  did  not  limit  the  number  of  points  that  could  be  used.  Therefore,  the  limit  was 
determined  by  round-off  errors,  and  a  double  precision  procedure  was  used. 

All  numerical  results  were  obtained  for  a  value  of  Poisson's  ratio  equal  to  0.  3. 
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Section  6 
RESULTS 


The  accuracy  of  numerical  results,  of  course,  depends  on  the  use  of  a  sufficient 
number  of  points  in  the  finite  difference  scheme.  On  the  other  hand,  the  computer 
time  increases  rapidly  with  an  increasing  number  of  points.  The  convergence  of 
the  method,  therefore,  was  explored  through  calculation  of  critical  loads  for  fixed 
shell  parameters  and  a  successively  increasing  number  of  points.  Some  of  the 
results  so  obtained  are  shown  in  Table  1.  By  use  of  these  exploratory  calculations 
it  was  possible  to  establish,  as  a  function  of  the  parameter  Z,  the  number  of  points 
needed  for  l/2  percent,  or  better,  accuracy  in  the  final  results. 

Computed  critical  values  of  the  axial  load  for  the  case  of  zero  lateral  pressure  are 
shown  in  Table  2.  For  comparison,  corresponding  results  were  also  obtained  by  use 
of  the  membrane  prebuckling  solution.  It  was  found  that  with  boundary  conditions 
corresponding  to  cases  S3  and  S4  (W  =  0  ,  N  -  0)  these  results  differ  very 
little  from  the  results  in  Table  2.  In  the  other  six  cases  the  critical  load  with 
membrane  prebuckling  solution,  within  the  parameter  range  considered,  is  equal 
to  or  insignificantly  higher  than  the  classical  buckling  load  for  simply  supported 
cylinders  ( N  =  1.0)  . 

It  appears  from  Table  2  that,  with  the  exception  of  very  short  shells  ( Z  -  50 )  ,  the 
critical  load  is  practically  independent  of  the  parameters  r/t  and  L/r  in  all  cases. 

In  contrast  to  expectations,  lower  values  of  the  critical  load  are  in  some  cases  found 
for  the  very  short  shells.  Therefore,  the  influence  of  the  shell  length  on  the  critical 
load  was  studied  in  more  detail.  The  critical  load  versus  the  parameter  L/r  ,  for 
a  cylinder  with  r/t  =  100  and  with  boundary  conditions  corresponding  to  case  C2, 
is  shown  in  Fig.  1.  Here  the  number  of  circumferential  waves  was  held  constant 
(n  =  8  )  .  For  long  shells  the  critical  load  is  independent  of  shell  length,  and  for  very 
short  shells  the  critical  load  is,  as  expected,  monotonically  increasing  with  decreasing 
shell  length.  In  the  intermediate  range  an  oscillatory  behavior  is  displayed.  It  may 
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Table  1  (Cont. ) 


L/(  2h) 

ncr 

Case  S3 

Case  S4 

120 

0.5054 

0.5038 

240 

.5025 

.5006 

360 

.5021 

.5004 

480 

.5019 

.5003 

600 

.5019 

.5003 

720 

.5019 

.5003 

(d)  r/t  =  1(T  ,  L/r  = 


Case  S3 

L/(  2h) 

ncr 

120 

0.679 

240 

.676 

360 

.6750 

480 

.6745 

600 

.6745 

720 

.6745 

)  r/t  =  10Z  ,  L/r  =  0.7  , 


L/(2h) 

ncr 

Case  Cl 

Case  C3 

60 

0.91015 

0.90997 

120 

.91019 

.91002 

180 

.91020 

.91003 

240 

.91021 

.91003 

300 

.91021 

.91003 

CRITICAL  VALUES  OF 
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be  seen  that  the  deep  minimum,  which  is  reached  just  before  the  critical  load  increases 
monotonously,  is  very  close  to  L/r  =  0.7. 

A  similar  variation  of  the  critical  load  with  L/r  was  found  for  the  cases  S3  and  S4. 

In  Fig.  2a  the  critical  load  is  shown  versus  L/r  for  case  S3.  It  is  seen  that  in  this 
case  the  general  behavior  does  not  change  when  the  influence  of  lateral  restraint  in  the 
prebuckling  analysis  is  neglected.  For  relatively  long  cylinders  the  curve  with  mem¬ 
brane  prebuckling  solution  is  slightly  below  the  rigorous  solution.  However,  for  very 
short  cylinders  this  difference  increases,  as  may  be  expected.  In  Fig.  2b  the  length 
dependence  is  shown  for  a  cylinder  with  boundary  conditions  corresponding  to  case  S4. 
In  both  these  cases  the  curves  are  valid  for  a  fixed  number  of  waves  (n  =  2)  . 

Critical  values  of  the  external  pressure  with  zero  axial  load  were  also  computed  for 
the  eight  different  sets  of  boundary  conditions.  In  this  case  the  results  agree  with 
those  reported  in  Ref.  7,  except  for  very  short  cylinders.  The  slight  discrepancy  for 
short  shells  was  observed  for  simply  supported  shells  also  in  Ref.  8. 

Critical  loads  under  combined  loading  are  shown  in  Figs.  3,  4,  5,  and  6.  In  cases  SI 
and  S2  (Fig.  3)  the  value  of  r/t  has  practically  no  influence  on  critical  combinations 
of  the  axial  stress  and  external  pressure  parameters,  within  the  range  of  geometrical 
parameters  under  consideration. 

The  analysis  with  the  membrane  prebuckling  solution  indicates  that  the  difference  in 
critical  loads  for  cases  S3  and  S4  is  negligible  for  all  values  of  the  pressure  param¬ 
eter.  However,  in  the  presence  of  internal  pressure  the  rigorous  solution  gives  dif¬ 
ferent  results.  This  is  shown  by  the  interaction  curves  in  Figs.  4  and  5.  Here  the 
number  of  waves  corresponding  to  minimum  critical  load  generally  is  two,  but,  for 
higher  values  of  the  pressure,  in  case  S3,  this  minimum  occasionally  occurs  at  a 
larger  number  of  waves. 

In  case  S3  the  curves  show  that  for  most  combinations  of  geometrical  parameters 
there  exists  a  range  of  the  pressure  parameter  within  which  three  solutions  are 
obtained.  Of  course,  when  the  axial  load  on  the  shell  is  increased  under  constant 
internal  pressure, only  the  lowest  of  these  solutions  is  meaningful. 

For  clamped  cylinders,  it  was  found  again  that  the  parameter  r/t  has  no  influence 
on  critical  combinations  of  the  axial  stress  and  pressure  parameters.  It  was  found 
also  that  tangential  restraint  at  the  edge  does  not  affect  the  critical  load.  Interaction 
curves  for  clamped  cylinders  are  shown  in  Fig.  6. 
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Fig.  4  Interaction  Curves  for  Case  S3  (cont. ) 


Fig.  5  Interaction  Curves  for  Case  S4 
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Fig.  6  Interaction  Curves  for  Clamped  Cylinders 
p  <  0,  102  <  r/t  <  104  (cont.) 


Section  7 
CONCLUSIONS 


The  results  obtained  here  are  in  agreement  with  the  results  of  the  example  considered 
by  Fischer,  and  also,  over  most  of  the  parameter  ranges,  in  reasonably  close  agree¬ 
ment  with  the  results  presented  by  Stein.  However,  for  very  thin  and  long  cylinders 
the  difference  between  present  results  and  those  presented  by  Stein  is  appreciable. 
Therefore,  a  computer  program  based  on  Eq.  (15)  of  Ref.  2  was  derived.  It  was  found 
that  if  a  sufficiently  fine  mesh  size  was  used,  the  numerical  results  were  in  good 
agreement  with  the  present  results  throughout  the  parameter  range. 

The  results  for  simply  supported  cylinders  with  the  edge  free  to  move  in  the  tangential 
direction  show  a  drastic  reduction  of  the  critical  load  in  comparison  with  the  classical 
value.  For  relatively  long,  nonpressurized  cylinders,  this  reduction  is  practically 
independent  of  whether  an  accurate  or  a  membrane  prebuckling  solution  is  used.  For 
pressurized  cylinders  and  for  very  short  cylinders,  the  rigorous  analysis  gives  higher 
values  of  the  critical  load. 

For  simply  supported  cylinders  in  the  case  of  zero  tangential  edge  displacement,  and 
for  clamped  cylinders,  use  of  the  rigorous  analysis  leads  to  a  reduction  of  the  critical 
load.  This  reduction  is  rather  moderate  in  comparison  with  the  discrepancy  between 
theory  and  tests.  It  appears  that  these  boundary  conditions  generally  are  applicable 
in  experiments,  and  therefore  the  edge  effects  studied  here  cannot  alone  explain  the 
discrepancy  between  theory  and  tests.  A  final  solution  of  this  problem  will  probably 
have  to  include  effects  of  initial  geometric  inaccuracies. 

The  results  for  cylinders  with  edges  free  in  the  tangential  direction,  first  obtained  by 
Stein,  are  somewhat  related  to  the  results  by  Nachbar  and  Hoff  in  Ref.  9.  In  Ref.  9 
the  cylinder  edge  was  considered  to  be  free  in  the  lateral  direction,  and  also  in  this 
case  the  critical  load  was  drastically  reduced.  Although  neither  of  these  conditions  is 
realized  in  experimental  analysis,  they  are  still  interesting  in  that  it  is  feasible  that 
in  practical  applications  elastic  restraint  at  the  edges  may  be  rather  weak. 
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